T97/036, SU-4240-660 



Effects of Self- Avoidance on the Tubular Phase of 
Anisotropic Membranes 



M. Bowick, 

Physics Department, Syracuse University, 
Syracuse NY 13244-1130, USA 

E. Guitter 

CEA, Service de Physique Theorique de Saclay, 
F-91191 Gif sur Yvette Cedex, France 



We study the tubular phase of self-avoiding anisotropic membranes. We discuss the 
renormalizability of the model Hamiltonian describing this phase and derive from a renor- 
malization group equation some general scaling relations for the exponents of the model. 
We show how particular choices of renormalization factors reproduce the Gaussian result, 
the Flory theory and the Gaussian Variational treatment of the problem. We then study 
the perturbative renormalization to one loop in the self-avoiding parameter using dimen- 
sional regularization and an e-expansion about the upper critical dimension, and determine 
the critical exponents to first order in e. 



email: bowick@npac.syr.edu, guitter@spht.saclay.cea.fr 



1. Introduction 



The statistical mechanics of isotropic tethered membranes has been extensively stud- 
ied [|l|J^. In a recent paper Radzihovsky and Toner (RT) showed that intrinsically 
anisotropic tethered membranes are surprisingly rich systems. In particular they exhibit 
an intermediate tubule phase between the crumpled and flat phases typical of isotropic 
tethered membranes [^H- The tubule phase is characterized by being extended in one 
direction and crumpled in the other. Furthermore any degree of anisotropy is expected to 
be relevant, so such systems could be widespread in nature and very important. It is not 
hard to imagine many situations in which the polymerization of a fluid membrane occurs 
anisotropically. 

Recently the existence of this tubule phase for physical anisotropic membranes has 
been confirmed by large-scale Monte Carlo simulations |T^] and the crumpled-to-tubule 
and tubule-to-flat phase transitions both observed. In the case of self-avoiding tethered 
membranes, current numerical evidence suggests that the crumpled phase is destroyed in 



physical dimensions This enhances the possible signiflcance of an ordered tubule 

phase for self-avoiding anisotropic physical membranes — the only transition left in this 
case may be the tubule-to-flat transition. 

In this paper we study the effects of self-avoidance in the tubule model of a self-avoiding 
tubule, previously introduced and analyzed by RT H. This model may be considered as 



the analog of the Edwards model of self-avoiding membranes [[r2| - [L4| , appropriately adapted 



to the tubular geometry, with bending rigidity in the extended direction of the tubule and 
self-avoidance in its crumpled direction. 

In Section 2 we use a renormalization group equation to reproduce some of the critical 
exponent scaling relations of RT and derive some new ones. These relations hold provided 
that the bending energy term is not renormalized, and imply that there is only one in- 
dependent exponent in the model. Special cases of this treatment reproduce the trivial 
Gaussian model as well as the Flory theory and the Gaussian Variational approximation 
results of RT. 

In Section 3 we establish the perturbative renormalizability of the model and prove 
that the bending energy term is indeed not renormalized. 

In Section 4 we calculate the critical exponents to flrst order in an e-expansion about 
the upper critical dimension for the relevance of self-avoidance. We use the techniques of 
dimensional regularization and the Multi-local Operator Product Expansion (MOPE) of 
reference [15|. We give the corresponding predictions of all relevant critical exponents for 
the case of a physical membrane in the tubule phase. 



2. Scaling Relations 

We start by reviewing RT's model and scaling results for the tubular phase of self- 
avoiding anisotropic membranes [^. We consider the generalized case of iI>-dimensional 
objects [D = 2 corresponding to membranes) with one stiff direction y and [D — 1) soft 
directions x± (see Fig. l-(a)). In the tubular phase, such an object will be extended in the 
y-direction and crumpled in the transverse direction. Using a Monge-like representation, 
the point with coordinates {x±, y) in the membrane will occupy a position {h{x±, y), y) in 




Fig.l: (a) An anisotropic membrane with a stiff y direction and a soft 
direction; (b) after embedding, the membrane forms a tubule, extended in 
the stiff y direction and crumpled in the soft direction. 

the (i-dimensional embedding space, with h a {d — 1) -dimensional vector field perpendic- 
ular to the y direction. Adapting the Edwards model for self-avoiding membranes to the 
geometry of the tubular phase RT obtained the Hamiltonian ^ 

^ = \j d^'-^x^dy (c>iK(xi,y))'| 



+ ^ J d^-^x±d^-^x'^dy 6^'^-^^ (h{x±,y) -h{x'^,y)^ . 



(2.1) 



The first two terms describe the elastic properties of the membrane in the absence of self- 
avoidance, and represent a bending energy term in the extended stiff y direction and an 
effective entropically generated elastic term in the crumpled direction. The third term is a 
two-body contact interaction with excluded volume (or self-avoiding) parameter b. Due to 
the extended nature of the tubule in the y direction, the self-avoiding interaction involves 
only points which have the same y coordinate along the membrane [Q. 

The engineering dimensions of the fields and coordinates are [y] = 1, [x^] = 2 and 

Co^[h] = ^-D. (2.2) 



This implies [b] = — e with 



SD---(--D]d. (2.3) 
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We consider the model for 3/2 < D < 5/2 only, where the bare roughness exponent Co 
( p.2|) satisfies < Co < 1- Setting e = fixes the upper critical dimension for the relevance 
of the self-avoiding interaction to be 0] 

with, in particular, dud!^) = 11. 

In Section 3 we will show that the Hamiltonian ( |2.1| ) renormalizes onto itself. In other 
words, one can find renormalization factors Z, Z± and Zf, such that the renormalized 
theory, using the renormalized Hamiltonian 



+ j d^-'xId^-'xfdyS^"-'^ {h^ix1,y)-h''{xf,y) 



(2.5) 



gives finite results at e = when expressed in terms of the renormalized self-avoiding 
parameter b^. We will moreover show that the bending energy term is not renormalized, 
that is 

Z = 1. (2.6) 

This assumption is crucial in the derivation of the scaling laws below. The bare ( |2.1| ) and 
renormalized ( |2.5| ) Hamiltonians can be made identical by appropriate rescalings of the 

height field h^, transverse coordinate x^ and self-avoiding coupling in the following 
way: 

h^{xU) = zt''^/'h{x^,y) 

x1 = Z]/^x^ (2.7) 

b^ = 6^-^Z-i4^-^)("+^)/^ 

Consider now the height fiuctuations in the bare model as determined by the correlation 
function 

Gix±. y) = -^^—^^ ( {kx^.y) - K(0, 0))'^ (2.8) 
From (12.71) the renormalized version of this correlation function satisfies 



G«(x^, y) ^ -^7^^ ( {h^'ixl, y) - /^^(O, 0))'\ = 4^-^)/'G(a:^, y). (2.9) 



R 



[z^f ^^^'^Gii{x^,y)^ = 0, where the derivative is taken at fixed bare 



Writing ^ 

parameter 6, we get the renormalization group equation 



where ^ = A* ^ 



logZ±. We suppose here that an infrared stable fixed point is reached, 



describing the large scale properties of the membrane. The equation ( p.lO|) holds precisely 
at this fixed point. On the other hand simple scaling gives the homogeneity equation 

fi^G^ - y^G^ - 2x^^G^ + (5 - 2D)G^ = 0. (2.11) 
on oy ox± 

We thus get the fixed point renormalization group equation 

^|_G'« + - 2CG« = 0, (2.12) 

oy z ox± 

where the anisotropy exponent z and the roughness exponent C are given by 

2 



4 + ^^ (2.13) 

C = Co + -^5 , 

with the bare roughness exponent Co given by (|2.2| ). Eq. (|2.12| ) implies the scaling 

Gnixf, y) ~ y^^F,{y/{xlY) ~ {xlf^F^{y/{xf)^), (2.14) 

where the size exponent v and roughness exponent C are related by = C^;. Eliminating 5 
in ( p.l3|) , we thus find the very general scaling relations 



, 3 l-D 

Iz il-D) ^2-1^) 



Rewriting Eq. ( |2.14| ) in momentum space and using the derived scaling relations one 



finds that the inverse of the height field propagator G~^{q-,p±) scales as G~^{q-,Pi) 

f{Q/p±^)- Thus the anomalous dimension, rj, for the bending rigidity vanishes, as re- 
quired by the non- renormalization theorem Z = 1. Similarly it is simple to show that 
G~^{q,p±) = p±'^~^^^g{q/p±^), with r]± = 4z — 2. Since the size exponent u must exceed 
its phantom value Co/2, one has z > 1/2 and therefore rj± > 0. 

From the above scaling relations, we end up with only one independent exponent in 
the theory, depending on the precise value of S. This value, and the subsequent predictions 
for all exponents, may be fixed by imposing one more constraint on the renormalization 
factors of our model Hamiltonian. At this stage, this extra imposed constraint is totally 
arbitrary and different constraints lead to different values of the exponents. It is interest- 
ing nevertheless to explore limiting cases where scaling is dominated by one component of 
the Hamiltonian only, either the elastic term or the self-avoiding interaction. The corre- 
sponding limiting values of the exponents indeed define the range of values in which the 
exact exponents are expected to fall. One can fix the scaling from the elastic terms only by 



assuming the absence of renormalization for the h field, i.e. by imposing Z± = 1, yielding 
(5 = 0. One then recovers the bare values C = Co? ^ = 1/2, = Co/2 of the Gaussian theory 
without self-avoidance. On the other hand we can consider the strong coupling limit where 
scaling is fixed by the self-avoiding term only. This yields the Gaussian Variational result, 
as discussed in Section 2.2 below and also treated in A third, intermediate, estimate 
of the exponents is the Flory result, obtained by balancing the elastic and self-avoiding 
contributions in the Hamiltonian, as discussed in Section 2.1 and in 0. Notice finally that 
these different estimates become exact and identical on the e = line, and can be used as 
starting points for a systematic expansion in the [D, d) plane around this line. This idea 
has been used in for the self-avoiding isotropic membrane Edwards model. In Section 
4, we will calculate the correction to the Gaussian, Flory and Variational estimates of the 
roughness exponent at first order in e and for fixed D = 2. 



2.1. Flory Theory 

In Flory theory one assumes that elastic energies are comparable to self-avoiding 
energies. If this is to remain true under renormalization one should require that both 
terms renormalize in the same way viz. Z± = Zi,. Given this assumption one finds from 
if7\) that 

6« = 6^-4^-^)('^+^)/'-\ (2.16) 
In this case the fixed point condition directly determines S in terms of e to be 



6f^ 



-4e 



{4+{D-l){d + 3)} 

with e as in (|2.3| ) . Using (|2.13| ) and ( |2.15| ), the size exponent u is then found to be 

(^ + 1) 



(2.17) 



{d+1] 



(2.18) 



which coincides with the Flory prediction found in [Q. This is nothing but the usual 
Flory result for a {D — 1) -dimensional self-avoiding object in a {d — l)-dimensional em- 
bedding space, and corresponds to treating the different transverse slices of the tubule as 
independent . The other exponents are likewise determined in this approximation to be 



Zf 



4 + {D-l){d + 3) 

Hd+1) 

3(D + 1) 
4 + {D-l){d + 3) 



[2.19) 



The corresponding values for the physical tubule {D = 2 and d = 3) are dp = —8/5, 
i^F = 3/4, ZF = 5/6 and Cf = 9/10. 



2.2. Gaussian Variational Approximation 



A different approximation one can make is to assume that the self-avoiding term is not 
renormahzed viz. Zi, = 1. This is exactly the approximation which is made in a Gaussian 
Variational treatment of the problem, where the exact density functional is approximated 
by the best possible Gaussian weight for the field h, using a variational principle 0. In this 
case the field h is renormahzed, but the self-avoiding interaction term is not. Repeating 
the above analysis in the case — 1, one finds easily from the fixed point condition 

(2.20) 



(2.21) 



{D-l){d + 3) ■ 

The size exponent u in this approximation, first obtained by RT, is 

^ 7(1^-1) 
^""^ (3d -5) ' 

and the other exponents are likewise determined to be 

_ (£)-!)((/ + 3) 

f"^"^^ (2.22) 

^"""^ " (dT3) ■ 

The corresponding values for the physical tubule are S^ar = —2/3, u^ar = 7/4, z^ar = 3/2, 
and Cvar = 7/6. The unphysical nature of these values {ly and C cannot exceed 1) indicates 
that, in this approximation, the tubule phase is unstable. For D = 2, in fact, one sees 
from ( p.21|) that the tubule phase in unstable below d = 4. It is known, however, that the 
Gaussian Variational method is a strong coupling method which usually over-estimates the 
size exponent. 



3. Renormalizability 

We now turn to the issue of the perturbative renormalizability of the theory for e > 0. 
We rely on the general formalism introduced in reference |jT^ for the treatment of non- 
local interactions. In the diagrams of the perturbative expansion in b we first identify 
the singular configurations of interacting points which contain possible divergences. We 
then use a short distance multi-local operator product expansion (MOPE) to analyze these 
singularities and show that they are proportional to the insertion of multi-local operators. 
A simple power counting argument allows us to extract from all singular configurations 
those which give rise to actual divergences. This, together with some symmetry arguments, 
singles out all the operators which require renormalization. From this analysis, we deduce 
that the Hamiltonian ( p.l| ) renormalizes onto itself, according to Eq. ( |2.5| ) and moreover, 
Z — 1, i.e. there is no renormalization of the bending energy term. Our analysis will be 
presented for L> = 2, but it could be easily extended to the range 3/2 < D < 5/2 where 
the roughness exponent Co = (5 — 2D)/2 satisfies < Co < 1- 



Let us concentrate on the partition function Z associated with the Hamiltonian 
at D = 2: 

= y V^{x,y)\ exp(-7^[K]) . (3.1) 
It can be expanded in powers of h according to 



Af=0 1=1 \i=l 



where Zq is the partition function of the non self-avoiding (6 = 0) theory and ((■■■))o 
denotes the corresponding Gaussian average 



((■■■))o = ^/ ^^[M^,2/)]exp(^-^ y dxdy | V (s.^x, y))' 



Each 8 function in (13.21) can be written as 



(3.3) 



5('^-i)(K(x.,|;.)-MxU.)) = I eMk^-y^)-H<,y.)) (3.4) 

and one is lead to evaluate the Gaussian average 



/ N \ AT 

/ JJe'^^-('^("-^^)-'^("^'^^)) \ =exp(^ - h- kj{Go{ 

\i=l I «,J = 1 



-Gq{x[- Xj.Vi-yj) (3.5) 
+6*0 (x- -x'j^Vi -Vj)} 



where Gq is the two-point functionll 

1 

'2(rf- 1 
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x\ ' exp — — — j- ] — -y erf 



(3.6) 



The term of order in the perturbative expansion ( p.2|) is therefore naturally represented 
by a diagram of "dipoles" of interacting points located at {xi,yi;x[,yi) with "charge" 
±/ci, as depicted in Fig. 2. Note that the two end points of a given dipole i are located 



^ Here erf(n) denotes the usual error function erf(n) = (2/y^) dtexp{—t^). 
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Fig. 2: The diagram of order in ( p.2| ) is made of dipoles. The two end 
points of a given dipole are located at the same position yi in the y direction 
but at different positions Xi and x'^ in the x direction. 

at the same position yi in the y direction but at different positions Xi and x'^ in the x 
direction. A singular configuration of interacting points is found when the quadratic form 

Q{{h}) = X^^i ■ kj[GQ{xi - Xj.yi - yj) - G'o(x- - Xj.yi - yj) - Go{xi - x^yi - yj) 



appearmg m 
function 



we obtain 



+ Go{Xi- x'j,yi-yj)} 
(3.7) 

is not positive definite. Using the integral representation of the two-point 



Go{x,y) 



dp dq 



dp dq e' ^P^+iy) - 1 
27r 27r q^ + p'^ 



2n 2n 



(3.9) 



The quadratic form Q is thus positive definite except for those configurations of end points 
{xi,x[,yi} for which one can find a set of charges {ki}, not all zero, satisfying: 



J2 k^ e^'^y^ (e^P"^ - e'P^'^ ) = V(p, q) 

i 

^ p{.x, y) = ^ki6{y - yi) {S{x - Xi) - 6{x - a;-)) = V(x, y) 



(3.10) 



This latter condition is the requirement that the charge density p{x, y) vanishes identically, 
while some of the charges ki remain non-zero. 

This is possible if some of the dipoles arrange to form a so-called "molecule" , i.e. attach 
their end points and assemble into a connected diagram with at least one loop, such as in 



X 




atoms 



^ y 

Fig. 3: A molecule with 2 loops made of a connected assembly of 4 dipoles. 
This molecule has 3 atoms located at different values of x but at the same 
value of y. 



Fig. 3. A set of end points in contact form what is called an "atom" and their common 
position is the position of the atom. A zero of Q is obtained by an appropriate choice of 
non-zero charges keeping all the atoms neutral, which is possible in the presence of a loop. 
Note that all the atoms of a molecule have the same position in the y direction. Note also 
that singularities coming from disconnected molecules can be treated separately and that 
dipoles in the molecule which do not belong to a loop (dead branches) do not contribute 
to the singularity and can be ignored. 

The construction above identifies the singular configurations of end points which give 
rise to possible divergences. Such a configuration is characterized by a set of M atoms 
labeled by p and with position (xp, y) with the same y coordinate for all the atoms. For 
each atom p, we denote by Ip the set of dipoles i which attach their first end point at the 
atom p (i.e. {xi.yi) = {xp,y)), and by Jp the set of dipoles j which attach their second 



end point at the atom p (i.e. 



[X, 



[X. 



,,!/)). The singularity can be analyzed by use of 



the general short distance Multi-local Operator Product Expansion (MOPE) introduced 
In practice, one can return to the operator level (left hand side of (|3.5| )) and write 



m 



15 



the contribution of the atom p in 



iki-h{xi,yi) 



n 



-ikj-h(x'j,yj) 



-ikj-h(x'.,yj) 



X 



iki-h{xi,yi) 



ikj-h{x'yj) . 



ielp 



j&Jp 



(3.11) 



i.e. separate in the right hand side of ( |3.5[ ) the propagators Gq which involve only points 
inside the atom p, and which reconstruct precisely the Gaussian average above, from those 
involving at least one end point not in the atom p, corresponding to a normal product 
prescription. This separation allows us to isolate the singularity in the factorized Gaussian 



average, while the normal product has a regular expansion in Xi — Xp {i & Ip), x'j — Xp 
{j E Jp) and yk - y (k E IpU Jp) 

. ^ik,M-P,y) (i+i (^^x^ _ a;^) k, . d^h{xp, y) + {y, - y) h ■ dyh{xp, y)^ 

-i 2^ (i^'j - ^p) kj ■ d^h{xp, y) + (y^ - y) kj ■ dyh{xp, y)J+...) : 

with kp = {J2iei ^» ~ '^je.J ^i) being the total charge of the atom p. The same treatment 
can be applied to all the atoms of the molecule, creating for each atom p an exponential 
factor e^'^p'^^i^P'y) ^ together with insertions of various and/or dy derivatives of the field 
h at the point {xp^y). As in the MOPE is obtained by performing the integration 



over the charges ki for the dipoles i forming the molecule. This expands the corresponding 
product of bi-local operators Hi (ji{xi,yi) — h{x[,yi)^ around the chosen singular 

configuration in terms of general M-body operators of the form 

„ M 

$(xi,...,xm)= d''-^hllAp(xp,y)vf'S^''-^^{h-h{xp,y)'^ , (3.13) 

p=i 



multiplied by singular coefficients (see for details). Here Ap{xp,y) denotes either the 
unity operator 1 or a local operator in the derivatives of the field h at point {xp,y) and 

is a short hand notation for na=i '^^^ above operators are multi-local in the x 

direction but local in the y direction. This is because all the atoms in the molecule have 
the same y-position. We will see two explicit examples of the MOPE in the next section 
where explicit one-loop calculations are presented. 

At this stage, let us mention the following important result concerning the case where 
one inserted operator involves dy derivatives only (such as {dyh)'^). Indeed, such a term 

comes from the expansion of some operator Q^^i-H^i^Vi) taken at Xi = Xp exactly (we 
suppose here that i E Ip rather than i E Jp). However, in contrast with the coordinate 
Xi (or resp. x'^) which appears only in the atom p, the coordinate yi appears in a second 
atom p' (such that i E Jp'), which is in general distinct from p. The expansion in yi ~ y 
can be done simultaneously on the operator Q^i<:i-h{xi,yi) ^i^ove for Xi = Xp and for the 
operator e"''^* '^^^*'^*) for x[ = Xp', in which case the operator to be expanded in yi — y 
is Q'^i-i^^'^p^yi)-''^^=^p'^yi)) , We thus get the important result that those operators with only 
partial derivatives in the y direction can be regrouped so that they involve the difference 
of the h field at two (in general different) points of the molecule. An example of such 
operator is the two-body operator: 

dy{h{xuy) - h{x2, y))y 5(^^-1) (h{xi, y) - h{x2, y)) (3.14) 

If the two end points of the dipole i happen to belong to the same atom p, then the operator 
to be expanded is 1, which means that this dipole cannot give rise to insertions of local 
operators with only dy derivatives. 



This latter remark has an important imphcation for the renormahzation of local op- 
erators, coming from the particular case of singular configurations where the molecule has 
only one atom. In this case, each dipole in the molecule falls automatically in the class 
just described of dipoles with their two end points in the same, unique, atom. We thus 
get the important result that local operators with only dy derivatives are not created by 
renormahzation. Such terms, when absent from the original Hamiltonian ( |2.1| ), never ap- 
pear, and the only such operator present in ( |2.1| ) (namely (dyh)'^) is not renormalized, 
that is 

Z = l. (3.15) 

Having identified the singular configurations and the corresponding general multi- 
local operators to which their singularities are proportional, it remains to identify those 
singularities which are not integrable and give rise to actual divergences. If the molecule 
is made of K dipoles, the operator which is expanded via the MOPE is the product of K 
5^'^~^^ factors, with canonical dimension —K(^o(d — 1) in units of y. The dimension of the 
multi- local operator $ in (|3.13|) is 



M 

{d - l)Co + - + - ' (3.16) 

p=i 

with the notation |m| = X]a=i "^a- "^^^ corresponding singular coefficient in the MOPE 
has thus dimension 

M 

-{K^l){d- l)Co - J](dim[Ap] - (|mp| + (c^ - 1)) Co) (3.17) 

p=i 

This coefficient has to be integrated over the 2K — M relative x coordinates and the K — 1 
relative y coordinates of the 2K end points of the dipoles approaching the positions of the 
M atoms. This gives a superficial degree of convergence for the corresponding integral 



2(2K - M) + (K - 1) + (M - (A' +!))(((- l)Co + $^(|nvlCo - diml^J) 

3(M - 2) + e (i^ - M + 1) + ^(ImpICo - dim[Ap]) 



with e = 5 — ((i — l)Co- Note that K — M + 1 is nothing but the number of loops in the 
molecule. A divergence is found whenever < 0. It is easy to check that all the local 
operators A but the unit operator have a strictly negative dimension in units of y, as a 
consequence of the relation Co<l- Ate = 0,Ci;<0 requires either M = 2, mi = m2 = 
and Ai = A2 = 1, which is nothing but the original contact interaction in ( |2.1| ), or M = 1, 
in which case $ is either the unity operator 1 or a local operator A{x, y) which moreover 
must satisfy dim[A] > —3. We already know from the previous discussion that A must 



contain at least one dj; derivative since terms with only dy derivatives are not created. Due 
to the X —X symmetry, the coefficient of a term with only one dj; vanishes and one thus 
needs at least two dj; derivatives. The term with largest dimension satisfying this criterion 
is the original elastic term in ( |2.1| ) (dxh)^ which already has dimension —3. It is thus, 
together with the unity operator, the only renormalized local operator in the theory. The 
renormalization of the unity operator is simply a shift in the free energy of the system. In 
particular, it disappears in the computation of average values of physical observables, and 
can be simply ignored. 

In conclusion, we have shown that the Hamiltonian 
Z = 1, as announced. 

For e > 0, the theory is super-renormalizable since uj in ( |3.18| ) increases with the 
number of loops of the molecule. 



2.1) renormalizes onto itself, with 



4. One-loop calculations 

Let us now present one-loop calculations, which give the corrections at first order in e 
for the critical exponents v and z. We use here dimensional regularization by considering 
the theory at e > and by calculating the renormalization factors and Z}, needed to 
make the theory finite for e = at one-loop order in 6^. We use a minimal subtraction 
scheme where we keep for the first order correction in Z^ and Z\j only the corresponding 
pole in e. 




(a) (b) 

Fig. 4: The two singular configurations leading to (a) a one loop renormaliza- 
tion of (dxh)"^ and (b) a one loop renormalization of 5*^^"^^ (jl{x, y) — h{x' , y)^. 

We obtain our results in the framework of the MOPE described above, which we use 
here in two simple cases: the 1 atom molecule made of a single dipole with its two end 
points approaching each other (see Fig. 4- (a)) and the 2-atom molecule made of two dipoles 
approaching each other (see Fig. 4-(b)). The first situation will give us a renormalization 
of the local operator {d^h)'^, and hence a first order correction to Z±. The second situation 
will give us a renormalization of the bi-local operator 5*^^"^^ y) — h{x' , y)j , and hence 

a first order correction to Zi,. 



To analyze the divergence for Fig. 4-(a), we use the Operator Product Expansion 
(OPE): 

^ik-{h{x,y)-h{x',y)) _ Go(x-x',0) 1^ ^ _ ^z^,^. . Q^^^^^ ^) . 

1 



where xq = {x + x')/2 and 



2 ,x - x'fko^kp : d^ho,{xo, y)d^h/3{xQ, y) : + ■ . 

(4.1) 



G'o(x-x',0) = ^\x-x'\^/^. (4.2) 

2v7r 



When integrated over k, it gives the MOPE: 

^('^-i) fMx,y) - Mx',y)) =-^\ ^ ^ X 1 



(47r)— '-(_G'o(x-x',0)) — 

: (^d^h{xo,y)^ : + ... 



4(-G'o(x-x',0))'^ 



(4.3) 



We use the above formula for the renormalized theory (^^) , expanded to first order in 
b^. As in (|3.2D , the bi-local 5 interaction comes with a factor and the singular- 

ity in ( [4.3| ) proportional to (dxh^)'^ will be cancelled by the corresponding counterterm, 
appearing with a factor —{Z± — l)/2, provided we choose: 

(^-1)-,.^/ (4.4) 



2 2 7|^.|<,.-,4 (4,r)= (|X|i/2)6 

where X — — x'^, and where we have used (|4.2|) and (i = 11 — 2 e. This leads to: 



6« 1 

Let us now analyze the divergence for Fig. 4-(b). We now use the OPE for the first 
atom 

^i[ki-h{xi,yi)+k2-hix2,y2)) _ g- fci •fc2Go (xi -X2 ,yi -y2 ) . gi(fci +fc2)-?i(a;o ,yo) ^ | . (4 6) 

where xq = (xi + X2)/2 and yo = (yi + y2)/2, and the similar OPE for the second atom 

^-i(^ki-h{x[,y-i_)+k2-hix2,y2)) _ g-^i •fc2Go(x'^ -Xj ,yi -^2) . +fc2)-/i(xo,yo) _j_ j . (4 7) 

where = {x'l +X2)/2. The MOPE is obtained by integrating over ki and k2- More 
precisely, we define k = ki + k2 and = {ki — /c2)/2, so that /ci • /c2 = — + 0(A;^). This 
latter (9(A;^) term can be set to zero if we are interested in the leading singularity, which 



is responsible for the divergence. Integrating over k reconstructs a bi-local 5 operator and 
we get the foUowing MOPE: 



{hix^, 2/1) - hix[, 2/1)) (/^(X2, y2) - h{^2^ y2] 



/A \ "^"1 d-l 

(47r) 2 _ - yi) - G'o(x';l ~ ^2' ?/i - 2/2)) ^ 

X 5(^-1) (Mxo,2/o)-Mx[„2/o))+.. 



(4i 



We are interested in the pole in e obtained when integrating the coefficient of the 5 term 
in the right hand side of ( ^^) over the relative coordinates xi — X2, x'l — x'2 and yi — y2- 
Defining Y = \yi — 1/2I + \xi — X2\^^'^ + \x'i — x''^!'^, u = \xi — X2\^^'^/\yi — 2/2! and 
V = — x^^/^/Iyi — 2/2I5 and using again the explicit formula ( p.6|) for Gq and d = 11 — 2e, 
we obtain a pole in e equal to: 



(4x)5 i„ ys- io i„ (/(«) + /(«>))' 



where /(m) = u exp(— l/4w^) + erf(l/2w). The integral over Y gives a pole "^/e. 

The integral over u and f is convergent and will be denoted by: 

poo pao 1 pea 

/-/ du dv-— ——-, = - da{F{a))\ (4.10) 
Jo Jo [.fM+f(v)) Jo 



where F{a) = a? duue~°' The function F{a) satisfies F{a) 1 and F{a) 
exp(— a^Tr/S). The integral / can be estimated numerically to / = 0.068373636(1). 

Applying as before the MOPE of Eq. ( [4.8|) to the renormalized theory, now expanded 
as in ( p.2| ) to second order in b^, the N = 2 diagram gives two d interactions with a factor 
{b^/j,'^)^ /8, leading to a divergence equal to 

2xMj,^l:z (4.11) 

8 TT 2 e 

with a factor of 2 coming from the two ways of assembling the two dipoles of the diagram 
into a one-loop molecule. This divergence will be cancelled by the 5 interaction counterterm 
in the renormalized Hamiltonian, which comes in the expansion with a factor — (Z^ — 
1)6^/^72, provided 

b^ I 

Z, = l + -^-. (4.12) 

2n2 e 



Using (|4.5|) and ( 4.12|) , we relate the bare and renormalized coupling constants as in 



for i:) = 2 and d = 11 - 2e: 



leading, after differentiation with respect to n at fixed b to the one-loop Wilson function: 







/ 7-e 1 
2^ ^ 2 16^ 



(6^)2 + C ((6^)3). (4.14) 



We thus obtain an infrared stable fixed point at 

e 



'R* 



l_ ± 

o S ^ 2 167r2 

ZTT 2 



(4.15) 



This fixes the value of the anomalous dimension 6 through 



log = logZ, = (-.(,« + O (((,«)^) ) (jg^i + O ((.«)) 



167r2 



and 



5 = 5(6 



^ I _ 



7 + 0{e' 



(4.16) 
(4.17) 



Numerical values for the exponents at D = 2 and d = 3 are obtained by setting e = 4 in 
the above formula, giving 

■1.050 (4.18) 



and thus the estimates 



z = 0.678 
u = 0.517 
C = 0.762 



(4.19) 



To understand the values we obtain for these exponents more clearly notice that the factor 
7/2 in the denominator of ( [4.17| ) is actually the factor {d + 3)/4 = (7/2) — (e/2), appearing 
in the exponent of Z± in ( |2.7D , to first order in e. It is therefore legitimate, at first order, 
to replace this factor 7/2 by the factor 3/2 obtained by setting d — 3 directly. This in 
practice amounts to making a partial two-loop correction. This leads to a new estimate 
5 = —2.212 and ( = 1.053, well above the original estimate ( |4.19| ), and actually even 
unphysical since larger than one. We see here that, due to the large value of e = 4 at the 
physical dimension d = 3, the first order estimates ( [4. 19] ) are not robust with respect to 
second order corrections and cannot be reliable. 

It is also interesting to develop alternative expressions for the roughness exponent 



as was done for the isotropic membrane Edwards model in |]T^. Indeed, the above estimate 
of C relies on the expression ( |2.13|) expressing the deviation of ( from its Gaussian value 



Co at e = 0, in terms of the anomalous exponent d = fx 



log Z± , which we estimated 



to first order in e in ( [4.17| ). Using the relation ( |2.7| ) between the bare and renormalized 



coupling constants, however, we can write, at the fixed point, the two foUowing equivalent 
definitions of 6: 



5 = 5f 



d 



{4+{D-l){d + 3)}'^ 



1 

log — 





d 



""^^ (D-l)((i + 3)^ d^ 
leading directly to the two identities 



(4.20) 



logZfc 



C = Cf + 



D-1 



d 



{4 + {D-l){d + 3)r 
logZfo . 



1 

^1. 



1 d 



(4.21) 



These relations express the deviation of C from its Flory value and its Variational value 
respectively. As we did for 5 in ( [4.16 ) and ( [1.17| ), we can obtain for D = 2 and d = 11 — 2e 
the estimates to first order in e: 



d 

dfx 



7 1 - ^ 

log|^ = ^e + 0(e^) 



8/ 



(4.22) 



log ^6 



8/ 



+ 



ye + 0(e 



One can easily check that the two equations in ( |4.21[ ) give exactly the same estimate as 
before for at first order in e, provided that the quantities Cf, Cvar and the different factors 
appearing in (|4.21 ), which involve d, are themselves expanded to first order in e. 

On the other hand, one could also decide not to expand any of these factors and 
impose d = ?> directly. If one moreover restores the factor (li + 3)/4 instead of 7/2, as 
discussed above, all the various expressions reproduce the unphysical estimate C = 1.053. 
If only some of the terms are expanded in e, we get lower values of C,. We thus expect 
that the original estimate C, = 0.762, obtained by expanding all terms at first order in e is 
actually a lower bound on the exact value of C,. 



5. Conclusions 

In this paper we have studied, within the e expansion, the effects of self-avoidance 
in the tubule model introduced by RT, going beyond their variational and Flory theory 
treatments of self-avoidance. We first show that the model is renormalizable and, fur- 
thermore, that the bending energy term is not renormalized. We then derive very general 
scaling relations for the critical exponents of the model at an infrared stable fixed point. 
These relations imply there is only one independent exponent. For special choices of the 
renormalization factors we are able to reproduce three different limits of the model viz. 
the trivial Gaussian model, the Flory approximation and the Gaussian variational approx- 
imation 0. This shows the power of this approach. We then treat the fluctuations of the 



model to one-loop in the self-avoiding parameter in an e expansion about the upper critical 
dimension. This yields predictions for all the critical exponents to first order in e. 

One should notice that our results have been obtained for an infinitely large mem- 
brane. For a finite membrane with extension Ly in the y-direction and L±_ in the transverse 
direction, finite size scaling laws can be derived in the above renormalization group frame- 
work [0. Due to the anisotropic nature of the tubular phase, however, there are many 
different scaling regimes, depending in particular on the relative scaling of Ly and L^. 

Finally, let us stress that the above analysis of renormalizability does not depend 
on the precise form of the Gaussian elastic term in the y-direction. One could imagine 
replacing the bending energy term by a tension term (dyh)^, describing for instance a 
tubule under longitudinal tension. The theory would then also be renormalizable in an 
e = 2D — 1 — {d — l){2 — D) /2 expansion, with again no renormalization of this tension term 
and only one independent exponent in the theory. In this case, however, the calculation 
cannot be performed at D = 2 directly where the upper critical dimension is infinite. 
As for self-avoiding isotropic membranes, a complete study of the problem for D < 2 is 
required. 

After this paper was completed we were informed by RT that the Hamiltonian Eq. ( |2.1J ) 
is not sufficient for a complete description of polymerized tubules in d=3. RT argued that 
a more involved Hamiltonian, including the anharmonic elastic terms of RT Eq. (5), in 
addition to the self-avoiding interaction, is needed. Since the present paper is rigorous 
and self- consistent we feel that it nevertheless makes a vital contribution to our present 
understanding of tubules. The analysis of the fuller model suggested by the remarks above 
presents a very definite challenge — to our knowledge there does not exist in the literature 
any proper renormalization group treatment of a theory with both nonlinear elasticity and 
two-body self-avoidance. 
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